Photoluminescence investigations of 2D hole Landau levels 
in p-type single Al^Gai^As/GaAs heterostructures 

M. KubisaQ and L. Bryja, K. Ryczko, J. Misiewicz 
Institute of Physics, Wroclaw University of Technology, 
Wybrzeze Wyspianskiego 27, 50-370 Wroclaw, Poland 

C. Bardot, M. Potemski 
Grenoble High Magnetic Field Laboratory MPI-FKF/CNRS, 
25 Avenue des Martyrs, F-38042 Grenoble Cedex 9, France 

G. Ortner, M. Bayer, A. Forchel 
Technische Physik, Universitdt Wiirzburg, Am Hubland, 97074 Wurzburg, Germany 

C. B. S0rensen 
The Niels Bohr Institute, University of Copenhagen, 
Universitetsparken 5, DK-2100 Copenhagen, Denmark 
(Dated: February 1, 2008) 

We study the energy structure of two-dimensional holes in p-type single Ali-^Ga^As/GaAs het- 
erojunctions under a perpendicular magnetic field. Photoluminescence measurments with low den- 
sities of excitation power reveal rich spectra containing both free and bound-carrier transitions. The 
experimental results are compared with energies of valence-subband Landau levels calculated using 
a new numerical procedure and a good agreement is achieved. Additional lines observed in the 
energy range of free-carrier recombinations are attributed to excitonic transitions. We also consider 
the role of many-body effects in photoluminescence spectra. 
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I. INTRODUCTION 

Valence band subbands in doped quantum wells and 
heterojunctions have complicated structure. Due to 
strong electric fields resulting from both band-edge dis- 
continuities and irregular charge distribution the heavy 
and light-hole states are mixed. This results in a strong 
nonparabolicity and anisotropy of energy levels and their 
nonlinear behavior under an external magnetic field. 
Hole subbands in doped heterostructures were studied 
experimentally, however much less intensively than con- 
duction band states. In magnetotransport measurements 
0, § H g I effective masses of holes were determined 
from Shubnikov-De Haas oscillations. It was shown jj], ||] 
that in single heterojunctions the spin degeneracy of va- 
lence band subbands is lifted as a result of the asymme- 
try of confining potential and the large spin-orbit cou- 
pling. The complex arrangement of hole Landau lev- 
els was investigated in cyclotron resonance experiments 
10' |(], @, [§j [§ • However the photoluminescence (PL) spec- 
troscopy in a magnetic field normal to the layers, very ef- 
fective in providing information on the conduction band 
electrons, was relatively seldom used for p-doped struc- 
tures |LC| [0], [L2[ [L3|, |l|, [15| |l§. Due to higher effec- 
tive masses both the cyclotron energy and mobility are 
smaller in the valence band by nearly an order of mag- 



nitude. As a result the ratio of Landau level broaden- 
ing (increasing with diminishing mobility) to their dis- 
tance becomes substantially larger, reducing the resolu- 
tion of measurements. Usually the studied luminescence 
resulted from band-to-band transitions. Volkov et al. 
investigated radiative recombination of two-dimensional 
(2D) holes with donor-bound electrons in specially de- 
signed heterojunctions where a monolayer of donors was 
introduced at some distance from the interface. In re- 
cent publications JlT], observations of subtle many- 
body effects like positively charged excitons or shake-up 
processes were reported in p-type quantum wells. To re- 
duce the broadening of PL features in these experiments 
free holes were created not by means of a modulation 
doping but using optical excitation. 

Calculations of hole Landau levels in p-doped het- 
erostructures are quite complex. Similarly as in the 
case of conduction-band subbands one has to solve self- 
consistently Schrodinger and Poisson's equations. But 
valence-band states are described by the 4x4 Luttinger 
Hamiltonian |2(| which considerably complicates the 

_ , |2l M ~ 



computations. In early attempts 21 
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hole eigenstates were expanded in a limited set of ba- 
sis functions and calculated by matrix diagonalization. 
The precision of such methods strongly depends on the 
choice of the basic set and is known to fail for larger 
magnetic fields unless the number of basis functions is 
increased [^l] . Bangert and Landwehr (2^] examined the 
importance of cubic anisotropy terms commonly ignored 
in calculations. Theory of exchange-correlation effects in 
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the valence band was developed by Bobbert et al. ptj ] 
who proposed a simple method of including many-body 
interactions in subband calculations. In Refs jl5| and 
p7| [2^1 improved numerical algorithms for solving the 
Luttinger Hamiltonian were elaborated. The structure 
of hole Landau levels calculated by different authors was 
compared with experimental results, however the agree- 
ment was rather poor. 

In this paper we investigate low-temperature photo- 
luminescence from p-type single AL^Ga^As/GaAs het- 
erostructures in a magnetic field perpendicular to the lay- 
ers. As it was first demonstrated by Yuan et al. |2j], Q 
the PL spectrum of a modulation-doped heterostructure 
is dominated by a broad line, named H-band, located 
in the energy range of bulk excitons and donor-acceptor 
pairs recombination. The characteristic feature of this 
line, observed both in donor and acceptor-doped struc- 
tures, is a strong blue shift with increasing excitation 
power. A tilted magnetic field experiment [pd[ proved 
that 2D carriers are involved in the H-band emission. 
Time resolved PL measurements [[H] displayed extremely 
long life times which strongly depended on the line po- 
sition. This was explained as a result of spatial sepa- 
ration (by the interface electric field) of carriers partic- 
ipating in the recombination process. Along with these 
observations the H-band emission was interpreted as a 
recombination of 2D particles confined near the junction 
(electrons or holes according to the doping) with pho- 
tocreated carriers of the opposite charge removed by the 
electric field far from the interface. In strongly doped 
heterostructures with high densities of 2D carriers the 
large interface electric field results in appreciable spatial 
separation of recombining electrons and holes. As a con- 
sequence the H-band emission has a free-carrier character 
which reveals in linear changes of the line position with 
a perpendicular magnetic field |ll], |l6[ [52[ In less 
doped structures an excitonic nature of the H-band line 
was observed as a diamagnetic shift of the transition en- 
ergy jl6| [m| . This clearly results from reduced separation 
of recombining particles which enhances their Coulomb 
binding. 

Free-carrier recombination observed in our samples 
with high densities of 2D holes enabled us to derive the 
structure of valence-subband Landau levels. In PL mea- 
surements we applied very low densities of excitation 
power (< 0.16 mW/cm 2 ), which increased the resolution 
and allowed us to discern the fine structure of H-band 
emission. In addition to free-carrier lines we observed a 
recombination of donor-bound electrons with 2D holes. 
The sample parameters as well as details of the experi- 
mental methods are presented in Sec. II. To analyze the 
results we calculated the energies and wave functions of 
2D hole Landau levels using a new exact numerical pro- 
cedure that will be described in Sec. III. We assumed 
that the non-doped GaAs layer is of p-type and the 2D 
holes form an accumulation layer at the interface. Usu- 
ally (see for example Ref. p5|) the background doping 
of pure MBE-grown GaAs is of p-type in the low 10 14 



cm' 3 range. This results in weaker electric field near the 
interface and consequently in somewhat lower 2D hole 
subband energies as compared with the previous calcu- 
lations, where the n-type residual GaAs doping was as- 
sumed. We also evaluated the matrix elements for optical 
transitions between valence and conduction-band Lan- 
dau levels. The theoretical results are compared with 
measured energies of PL transitions in Sec. IV and a 
very good agreement is established. Besides the free and 
donor-related transitions new lines are discerned in the 
spectra at higher magnetic fields, possibly of excitonic 
nature. We discuss the influence of exchange-correlation 
effects on hole Landau levels and observed recombination 
processes. Section V will summarize our conclusions. 



II. EXPERIMENTAL SYSTEMS AND SAMPLE 
PARAMETERS 

We investigated two single p-type GaAs/Gao.5Al . 5 As 
heterostructures with hole concentrations pi = 7.6xlO n 
cm' 2 and P2 = 9.8 xlO 11 cm' 2 . They were grown by 
molecular beam epitaxy on (001) semi-isolating GaAs 
substrates in the following sequence: 1 mm undoped 
GaAs, 7 nm undoped Gao.5Alo.5As spacer, 50 nm 
Gao.5Alo.5As layer doped with Be (5x 10 17 cm' 3 for sam- 
ple 1 and lxlO 18 cm -3 for sample 2) and 5 nm GaAs 
cap layer also doped with Be (lxlO 18 cm -3 for sample 
1 and 2xl0 18 cm -3 for sample 2). Transport properties 
of both structures in helium temperatures were reported 
previously P, |5[ |36[. 

Photolumincscence measurements were performed us- 
ing two set-ups. For preliminary studies we used an opti- 
cal split-coil cryomagnet providing magnetic fields up to 
8 T. The samples were excited in the Faraday configura- 
tion with Ar + laser (wavelength 488 nm) . Luminescence 
spectra were analyzed in 1 m single grating monochroma- 
tor with Nitrogen-cooled CCD camera used as a detector. 
The second set-up allowed extending PL measurements 
up to the magnetic field of 14 T. An optical fiber system 
was used with the quarter-wave plate and a linear polar- 
izer placed together with the sample in liquid helium. In 
order to examine the circular (a + and a~) polarization 
of photoluminescence the direction of magnetic field was 
changed. The spectra were analyzed in the same way as 
in the first set-up. 



III. THEORY 

To calculate the energy spectrum of 2D holes we de- 
veloped a technique that can easily be explained for a 
particle with parabolic dispersion. Suppose that we want 
to determine bound-state energies and wave functions of 
the Hamiltonian 
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in which both the potential V and effective mass m are 
position dependent only inside some interval Zi < z < 
Z 2 . By introducing points zq = Z\ < z\ < . . . < zjq = Z2 
we divide this interval into N parts and then approximate 
the actual functions V(z) and m(z) in each subinterval 
by appropriate constants (such as Vj = V(zj) and rrij — 
m(zj), respectively, for Zj—i < z < Zj). After that local 
solutions of the Schrodinger equation can easily be found: 

Fj{r) = exp(ikxrj_)fj(z) (zj-i < z < Zj), 



fj(z) — Aj exp(ikjz) + Bj exp(-ikjz), 



where 



2m j 



(E-Vj)-kl, 



(2) 



(3) 



E and k± = (k x ,k y ) are the particle energy and wave 
vector, respectively, and Aj, Ba are some unknown coef- 
ficients. Equations (||) and (g) are also relevant in the 
exterior regions z < zq (j — 0) and z > zn (j = N + 1). 

The matching conditions (i.e. the continuity of both 
the wave function and its derivative divided by the effec- 
tive mass) at points of division can be handy written in a 
matrix form. Introducting column vectors Uj 
we obtain 



B j-i( z j)Uj-i = Bj(zj)Uj, j = 1, 
where 



,N+1, (4) 



Bj(z) 



exp(ikjz) exp(—ikjz) 
£±exp(ikjz) -^exp(-ikjz) 



(5) 



Equations @) allow to relate the vectors Uo and Un+i in 
the exterior regions: 



A 
Bo 



T 



An+i 
Bn+i 



(6) 



where the transfer matrix T is given by the product 

N+l 

T =Il B iM*j-i)Bj(*3)- (7) 

3=1 

For a bound state the wave vectors A: and fc^+i must be 
imaginary as well as the coefficients Aq and Bn+i have 
to vanish. According to equation (0) this implies 



T n {E)=0. 



(8) 



To determine the energies of bound states we calculate 
numerically transfer matrix (^|) for different values of the 
variable E (and a fixed wave vector k±) searching for 
zeros of the diagonal element T\\. After that the wave 
functions can be evaluated by assuming Ajy+i = 1 and 
calculating the coefficients Aj , Bj in successive subinter- 
vals from equations (0) . The precision of results depends 



on the density of division (number N of subintervals) 
and can easily be controlled. If the potential V(z) is not 
constant for z — > ±00 a good approximation of accurate 
bound states can be obtained by choosing the outermost 
points of division [zq and Zn) suitably far from the clas- 
sical turning points. 

This method of calculations can be also applied to ma- 
trix Hamiltonians. Holes in the valence band of III-V 
compounds are described by the Luttinger matrix j20[ 



Hr = - 



P+Q L M 

L+ P-Q M 

M+ P~Q -L 

AI+ -L+ P + C 



(9) 



where 



P 



M 



2rriQ 
L = 



K 2 ,Q- 

-- -n/3 



2m^ K * +K «- 2K ^ 
{K x 



mo 



2mn 



K 



-iV, 



-Ky)K z , 

iVi — —K X K V 
mo 



(10) 



too is the free electron mass and 71, 72, 73 are Luttinger 
parameters. The eigenenergies of Hamiltonian (^) are 

E±{k) = -P(k)± 



\jQ 2 {k) +L(k)L+(k)+M(k)M+(k), (11) 

where E + and E_ relate to light and heavy holes, re- 
spectively. The wave functions can be written as F(r) = 
Fjk)exp(ikr), where F(k) is one of the following vectors 



Li(k) = 



Hx{k) 



r Ri 




- ■ 


L+ 


,L 2 (k) = 


M 


M+ 


-L 







. Ri . 


—L - 




- -M - 


R2 



,H 2 (k) = 




Ri 


-M+ 




L+ 



(12) 



Here R 1 = Q - P - E + , R 2 = Q + P + £L, vectors L x , 
L 2 correspond to light-hole states and vectors H\ , H 2 to 
heavy-hole ones. 

In heterostructures the Luttinger parameters have dif- 
ferent values in subsequent layers. Thus both compo- 
nents of the Hamiltonian Hl + V(z) describing hole sub- 
bands vary with the distance z. To find the subband en- 
ergies we proceed similarly as for Hamiltonian (|]J). The 
local solution of Schrodinger equation in the subinterval 
Zj-\ < z < Zj, where values of V, 71, 72 and 73 are set 
fixed, can be written in the form 



Fj(f) = exp(ik±rx)fj(z), 
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fj(z) = [CiH x {k H ) + C 2 H 2 {k H )\ exp(ik„z) + [C^L x {k L ) + C 4 L 2 (k L )} exp(ik L ; 



[C 5 Hi(-k H ) + C 6 H 2 (-k H )} exp{-ik H z) + [C 7 Li(-fc £ ) + C 8 L 2 (-k L )] exp(-ik L z). 

I 



(13) 



Here the complex wave vectors fcj,, ku are roots of equa- cent subintervals can be related using the matching con- 



tions 
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E + (k ± ,k z )=E-V j , E_(k±,k z ) = E-V j , (14) 



respectively. We choose k^ and ku such that their imagi- 
nary parts are nonnegative. Function (13) contains eight 
unknown coefficients, which will be represented by the 
column vector Uj = [Ci; C 2 ; . . . ; Cg]. Vectors Uj in adja- where Dj is the matrix operator 



fj-l(Zj) = fj(Zj), 



(15) 



D 



(71 - J 2 )d/dz V3"/3(k x - ik y ) 
-V3l3(k x + iky) {ji+j 2 )d/dz 



with values of Luttinger parameters appropriate for the 
subinterval j. Once more we obtain equations (0), this 



J 






(71 +l 2 )d/dz 
a/373 (&x + ik y ) 



-VS^ikx - iky) 
(li-72)d/dz 



(16) 



time with 8x8 matrices Ba of the form 



B, 



Li(fci) exp{ikLZ) L 2 (kL)exp>(ikLz) iJi(fc#) exp(ifc#z) H 2 (ku) exp(i/cj^z) 
DjL\(kL) exp(i/ciz) DjL 2 (kj}) exp(i/c£z) DjHi(kn) exp(ifc#z) DjH 2 (kn) exp(i/c#z) 



L\(-k L ) cxp(-ik L z) L 2 (-k L )cxp(-ik L z) Hi(-k H ) cxp(-ik H z) H 2 (-k H ) exp(-iknz) 
DjLi(-kL) exp(— ifciz) DjL 2 {— ki,) exp(— ik^z) DjHi^—kji) exjp(—iknz) DjH 2 {—kH)sxp{—ikHz) 



r 



(17) 



To find an equivalent of condition (JSj) for bound-state 
energies we examine the structure of vectors C/o and Un+i 
in the exterior regions. They are related by the equa- 
tion Uq = TU n+1 with transfer matrix (0) evaluated by 
means of matrices (^). The complete wave function is 
evanescent as z — ► ±00 only if imaginary parts of both 
/cl and ku are nonzero (i.e. positive according to the 
definition) for z < zq and z > zn- Moreover the compo- 
nents C\, . . . , C4 of vector Uq must be zero and similarly 
C5 = C*6 = C7 = C$ = in the case of vector £/at + i as 
can be deduced from expression (13). We thus obtain the 
equation 



where 



Ti 







■o- 


c 2 







c 3 







.c 4 . 




.0. 



(18) 



Ti = 



Tu ■ ■ ■ T14 
T41 . . . T44 



is a 4x4 submatrix of transfer matrix T. 
solutions of equations ( fl8| ) exist only if 

det[Ti(£)] = 0, 



(19) 



Non trivial 



(20) 



which is the requested condition for energies of bound 
states. To calculate the wave functions we first deter- 
mine the components C\ , . . . , C4 of vector U n+i by solv- 
ing equation (|l8|). After that relations (||) are used to 
compute vectors Uj in succeeding subintervals. Numeri- 
cal calculations proceed similarly as in the case of scalar 
Hamiltonian (|l|) except that higher-order matrices (ITT 
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must be evaluated. 

To account for a magnetic field B parallel to z-axis 
operator K (O) is replaced by K — — iV + eA/hc and 
the k terms p8| are included in Hamiltonian (^). Here 
A is the vector potential, e is the electronic charge and 
c is the velocity of light. The resulting Hamiltonian can 
be written in terms of the creation and destruction op- 
erators a = \/hc/2eB(K x ± iK y ). We apply the ax- 
ial approximation |E4| by replacing operator M ( [To| ) by 
M = V3(fr 2 j/2m ){K x ~ iK y ) 2 , with 7 = ( 7a + 73 )/2. 
Then the column vector 



*i(z)^n-l(r±) 
F 2 (z)4> n (f ± ) 

F 3 (a0&H-l(*L) 

•F 4 (z)0 n+2 (7\L) 



(21) 



becomes an eigenvector. The harmonic oscillator func- 
tions 4> n satisfy a + <f) n — y/n+ l<j> n +i and a~(f> n — 
y/n<j> n _\. In expression ( ^l| ) the index n runs over 
the values n = —2,-1,0,1,... and the envelope func- 
tions Fi(z) are automatically zero for those components 
which have harmonic oscillator functions 4> n with n neg- 
ative. By substituting vector ( pl| ) to the Schrodinger 
equation H^F = EF we obtain a set of ordinary 
differential equations for the z-dependent components, 
which can be rewritten as some 4x4 matrix acting upon 
F\{z), . . . , F^(z) functions only. The solutions of this 
set can easily be obtained with a resemblance to equa- 
tions (jll]) and (|j~2|). The procedure of transfer matrix 
evaluation following from this point is analogous to that 
for the B = case. 

In a selectively doped single heterostructure each hole 
moves in a potential which can be written as a sum 



V(z) = V (z) + V I (z) + V s (z). 
Here the component 



V (z) = 



-V , z^Q, 
0, z > 0, 



(22) 



(23) 



describes the valence band offset at the junction (z = 0). 
The potential Vj(z) created by the immobile charge (ion- 
ized impurities) and the potential Vs{z) due to mobile 2D 
holes can be both calculated in the Hartree approxima- 
tion by solving Poisson's equation 



d 2 V 

dz 2 



4ne 

—p(z) 



(24) 



with the suitable charge density p[z). The dielectric con- 
stants e of GaAs and Gai-^Al^As are slightly different, 
but the resulting image potential was shown [ p9[ to have 
quite small effect on subband energies. Therefore we ne- 
glect it in the calculations. 

Our structures are doped with acceptors of density 
Nab m the Ga 1 _ 2 . Al^ As barrier, however a layer of thick- 
ness w (the spacer layer) is left undoped to improve 2D 
carrier mobilities. The holes removed from the acceptor 



levels form the 2D carrier gas with areal density Ns as 
well as the accumulation layer of width I a- We assume 
that all barrier acceptors are ionized for —lj < z < —w 
and similarly all GaAs donors for < z < I a- Then the 
resulting immobile charge density is 



—eNAB ( — h < z < —w) 
+eN D (0 < z < l A ) 



(25) 



where No is the concentration of residual (minority) 
donors. Residual acceptors with concentration N A (> 
No) are neutral in the accumulation layer. The condi- 
tion of charge conservation can be written as: 



NabIi 



Nc 



Nnl 



D>>A- 



(26) 



Calculating the density of mobile charge we assume that 
only the ground valence subband is occupied by holes. 
As we will see this assumption is fulfilled even in our 
highly doped structures. We also neglect a coupling of 
the hole motion in the z direction and the motion in the 
plane of junction, assuming that all 2D holes have the 
same envelope function f(z). Then the density of mobile 
charge is 



p s (z) = eN s f 2 (z). 



(27) 



Equation ( |24| ) with charge densities (g5|) and ( p7| ) can 
easily be solved and after some manipulations the results 
can be written as: 



Vi(z) 



2ne 2 



-N s z ~ 2{N acc - 
(2AW + N s )z 
N s z + 2N acc z(l 
^ N s z + N acc l A 



Ns)w [z < — w) 

( - w < z < 0) 

- z/2l A ) (0 < z < l A ) 
\z > Ia) 



Vs{z) 



2TT6 2 



-Ns\z- 



- 2z / f 2 ( Zl )d Zl -2 ( Zl - z)f(z 1 )dz 1 \. 
Jo Jo 

(28) 

where N acc — NoIa is the areal density of charged impu- 
rities in the accumulation layer. In formula (pq) we used 
condition (^6|) and neglected for z < — w the component 
of potential Vi (z) square in z (which has a small effect on 
subband energies since the hole envelope functions vanish 
rapidly in the barrier). The potential value at the end of 
accumulation layer: 



V(Ia) 



2ne 2 



■Nn 



Ane 2 



N s I z^ 2 {z x )dz x 




(29) 



defines the position of valence band edge far from the 
interface (in the flat-band region) and is related to the 
hole Fermi energy Ep. If we denote by Ea the position 
of acceptor levels outside the accumulation layer and by 
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TABLE I: The valence-band parametres used in the calcu- 
lation (after Ref. [^o|). Linear interpolation is applied for 

AlarGai-aAs. 





GaAs 


AlAs 


7i 


6.85 


3.45 


72 


2.10 


0.68 


73 


2.90 


1.29 


K 


1.2 


0.12 



ea the bulk acceptor binding energy, we can write (see 
Fig. [l]) V(Ia) — Ea+£a- Solving the neutrality condition 



N. 



D 



N A 



1 + 2 exp( 



E; 



with respect to Ea we obtain 



V(l A ) = E F + e A + k B Tln 



N A -N L 



2N 



D 



(30) 



Equations ([29|) and (|30j) allow us to calculate both the 
length I a and density N acc of the accumulation layer. 
The Fermi energy Ep is determined from the known con- 
centration of 2D holes. 

Since the potential Vs(z) produced by 2D carriers de- 
pends on the hole envelope function, the subband calcu- 
lations must be self-consistent. At the first step we com- 
pute the potential V(z) and the energy and wave function 
of the highest heavy-hole subband edge. After any itera- 
tion the Fermi energy Ep is determined and new values 
of I a and N acc are calculated. If the in-plane wave vector 
k± is equal to zero, all off-diagonal elements of Hamilto- 
nian (^ vanish. Therefore the subband edges can be cal- 
culated exactly by solving the Schrodinger equation with 
scalar Hamiltonian (|IJ) and effective masses 7710/(71 — 272) 
for heavy and 7710/(71+272) for light holes. The Luttinger 
parameters used in the calculation are given in Table |. 
The value of valence-band offset Vb = 219.2 meV is ob- 
tained as 35% of the difference in band gaps in two adja- 
cent layers. The concentration of residual donors Np> is 
taken to be 5xl0 14 cm -3 . The concentration of residual 
donors Np> is taken to be 5 x 10 14 cm -3 . 

Figure |l| shows the computed valence band profile and 
positions of the first three subband edges for the sample 
with hole density N$ = 7.6xlO n cm~ 2 . Our results can 
be directly compared with those obtained in Ref. ||, 
where hole levels in the same structure were calculated 
assuming the n-type residual GaAs doping. In our case 
the electric field in the well region is substantially (almost 
two times) weaker. As a consequence we obtain higher 
values of hole subband edges and smaller intersubband 
distances (HH1 = —38.1 meV and LH1 = —45.9 meV as 
compared with values HH1 w —42 meV and LH1 as —56 
meV reported in H). In spite of a large 2D hole density 
only the ground HH1 subband is occupied and the Fermi 
level is located about 3 meV on top of the next LH1 sub- 
band. 



z(nm) 
100 150 




FIG. 1: Calculated valence-band profile in the well layer of a 
Alo.5Gao.5As/GaAs heterostructure with the 2D-hole density 
of 7.6 x 10 11 cm -2 . HHn and LHn label the edges of heavy and 
light-hole subbands, respectively. The dotted line represents 
the acceptor level. 



Using the determined self-consistent potential V(z) we 
next solve the Schrodinger equation (Hl + V)F = EF 
and calculate the subband dispersions for k± 7^ 0. Re- 
sults are plotted on Fig. ||. The subbands are spin-split 
(as a consequence of asymmetry of the confining poten- 
tial) and strongly nonparabolic. A comparison with the 
results of Ref. || reveals that the difference in potential 
profile affects also the in-plane hole masses. For exam- 
ple a hole in the upper spin-split subband HH1 is in our 
model almost 30% lighter. 

Calculations including the magnetic field are performed 
using the same self-consistent potential as for B = 0. We 
therefore neglect any influence of the magnetic field on 
the band bending. The results will be discussed in next 
section. 

To examine measured photoluminescence spectra we 
need to derive selection rules for interband optical tran- 
sitions. The complete hole wave function for a given Lan- 
dau level n can be written as 

= Fl4> n -lU\ + F 2 (j> n U 2 + F 3 <j) n+ lU3 + F 4 0„ + 2M4. 

(31) 

where it; (I = 1, 4) is the periodic part of the Ith 
Bloch function. The basis for Hamiltonian (||) is |37| 



U3 



(X + iY) i -2Z1 



1 

71 



(X - iY) T +2Z i J , u 4 = -j={X - iY) i . 



(32) 



Similarly the total wave function of a conduction band 
electron has the form = F e (z)<fi n u e where u e = S f 
or u e = S I and the envelope function F e (z) represents 
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FIG. 2: Energy bands as a function of the wave vector k± per- 
pendicular to the growth direction. The sample parameters 
are the same as in Fig. hi 



TABLE II: Matrix elements of allowed interband transitions 
involving the lowest conduction-band Lnadau level for two 
different circular polarizations. 



polarization 


hole 


electron 


matrix 




Landau 


spin 


element 




level 








n = 1 


T 


P f F?(z)F e (z)dz 


n = 


! 


(P/V3) f F^(z)F e (z)dz 


a+ 


n = -1 


T 


(P/V3) [FS(z)F e (z)dz 


n = -2 


i 


P j Fi(z)F e {z)dz 



in 



B=10T 




H 


3 


o 


D3 
I 


H2 




o 


D2 A 

W 


I 




D-A 


D1 / I 
••/ \ 


H1 / E1 
\J 


A°X 

UTf 

AIJLL. 



1.49 1.50 1.51 1.52 

Energy (eV) 



an unbound motion. We calculate the matrix element 
(^>n \ep\ ^ e n i) for the transition in which the light with 
polarization e is absorbed and an electron is raised from a 
state 9% to 'F^, . Since among matrix elements (S \p v \ X), 
(S \p v \Y), (S \p v \Z) (v = x,y,z) merely 

(S \ Pv \X) = (S\p v \Y) = (S \ Pv | Z) = im P, (33) 

are different from zero (P being the Kane momentum 
matrix element), we obtain 

(*£ \ep\ F e {z)^ n ,S T) - -im p\e- f Ff^F^d 3 ? 



and 
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FIG. 3: Photoluminescence from sample with Ns ~ 7.6 x 10 11 
cm -2 measured at B — 10 T in different circular polarizations 

2 



(34) 



for the excitation power density of 0.16 mW/cm 
indicated by arrows are identified in the text. 



The lines 



lowest conduction-band Landau levels with n' = partic- 
ipate to the luminescence. An inspection of formulas ( [34] ) 
and (|35|) gives the selection rules that are summarized 
in Table O. We consider only optical transitions in the 
Faraday configuration. The transitions involving valence 
band Landau levels n — and n — — 1 are three times less 
intense than the others, i.e. they have light-hole charac- 
ter. One can expect that the strongest transition takes 
place between the hole level n — — 2 and the spin-down 
electron state in the a + polarization. All the components 
of hole wave function ( pl| ) different than F^(z) (the ac- 
tive component according to Table [n]) are then equal to 
zero. 



(^ h n \ep\F e {z)4> n ,S i) = -m P 



+E 1 



F^* n+2 F e ^, n ,d 3 f+i 



F*<P n F e 4> nl d 3 r 



F;<t>*F e <i> n ,d 3 r 



(35) 



where s = (e x ± ie y )/^/2. In our p-doped structures the 
number of photoexcited electrons is small and only the 



IV. RESULTS AND DISCUSSION 

Photoluminescence spectra from the sample with A^g — 
7.6 x 10 11 cm -2 measured at B = 10 T in both a + and a~ 
polarizations are shown in Fig. |^. At low excitation pow- 
ers most of the incident light is absorbed in the GaAs 
region near to the interface. This explains relatively 
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FIG. 4: Shift and splitting of selected photo-luminescence lines 
in a magnetic field applied perpendicularly to the interface. 



small intensities of the bulk lines named X, D°X and A°X 
and attributed to free, donor-bound and acceptor-bound 
GaAs excitons, respectively. Stronger peaks denoted as 
H and D have two-dimensional character. This appears 
from the lines behavior with ch ang ed excitation power 



examined in our previous paper 1 33 . We found that the 



relative intensities of X-lines increased as the laser power 
was raised and more of the incident light reached the hat- 
band GaAs region far from the interface. At the same 
time the X-line positions remained unchanged whereas 
those of both H and D peaks shifted towards higher en- 
ergies. This results from the generation of electron-hole 
pairs which screen the electric field inside the accumula- 
tion layer. In effect the band bending is reduced and the 
energy of electron near the junction increases. On the 
other hand recombination processes in the flat-band re- 
gion remain unaffected. The line named D-A arises from 
the donor-acceptor recombinations inside the accumula- 
tion layer. Both its position and intensity also depend 
on the laser power. 

We identified the H lines as resulting from the recombi- 
nation of photoexcited electrons with 2D holes confined 
close to the junction. The lines named D were recognized 
as originating from the recombination of 2D holes with 



electrons bound to residual donors |33|. Positions of all 
free-carrier transitions and of two donor lines with lowest 
energies as function of the magnetic field are plotted in 
Fig. |j (the lines denoted as E will be discussed later). 
They display different behavior: the D-lines exhibit dia- 
magnetic shift whereas the energies of free-electron tran- 
sitions change almost linearly with B . At low fields 
the donor peaks are red-shifted with respect to the H- 
lines by the value of 4 meV close to the bulk GaAs Ry- 
dberg energy and the difference increases with growing 
magnetic field. 

To attain a high resolution the PL measurements were 
performed using very low density of excitation power 
P = 0.16 mW/cm 2 . We observed that at excitation 
levels below 0.7 mW/cm 2 the positions of H-lines sat- 
urate and do not change with the power P. Presumably 
the number of photocreated carriers is so small that the 
band profiles remain unaltered. Under these conditions 
the donor-related emission is distinct even at zero mag- 
netic field but the H-lines emerge only in higher fields 
B > 2 T. With increasing excitation power the donor 
lines gradually disappear [[33| . This can be related to the 
known effect of saturation of a defect-related photolumi- 
nescence Q. 

In this paper we analyze only the free-carrier transi- 
tions in an attempt to reveal the structure of valence- 
band states. To compare experimental results with 
calculated hole Landau levels we subtracted the mea- 
sured transition energies from the energy E cv + 0.5fiw c ± 
0.5/iBff-B of a GaAs conduction band electron in the low- 
est Landau level. Here \ib is the Bohr magneton, g 
= —0.44 is the electron g-factor and the effective mass 
m e — 0.067m was used to calculate the cyclotron fre- 
quency lu c . For each transition the projection of elec- 
tron's spin was resolved according to Table II. The energy 
E cv was used as a fitting parameter and a good overall 
agreement was achieved for the value E cv = 1464 meV 
(which will be discussed later). 

The results are shown in Fig. |[ We plotted merely the 
positions of Landau levels with highest energies in the 
valence subbands HH1 and LH1. The n = —2 level ap- 
pears only in the subband HH1 while Landau levels with 
ri ^ — 1 can be found in both subbands and for n ^ +1 
there are two levels with the same number n in each sub- 
band (n = level emerges once in the subband HH1 and 
twice in LH1). At lower magnetic fields all the transitions 
displayed in Table II are observed. They successively dis- 
appear with increasing field as the involved hole Landau 
levels cross the Fermi level. The line H5 related to a hole 
from the lower n = 1 Landau level (in the HH1 subband) 
vanishes at B w 3.5 T. This magnetic field value corre- 
sponds to the filling factor v = 10 at which the hole level 
becomes empty. 

The next line H4 observed in the o~ polarization in- 
volves a n = hole. It disappears somewhat below B = 6 
T and simultaneously a new line, denoted as El, emerges 
in the at slightly lower energy. This can be seen in detail 
in Fig. ^ where we show the PL spectra recorded at fields 
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FIG. 5: Hole Landau levels in the sample with Ns = 7.6 x 10 11 cm -2 calculated (lines) and measured experimentally (symbols). 
In the ground subband only the levels are shown which influence observed transitions and optically active hole states are 
distinguished by thick lines. Squares and triangles represent emission lines detected in the a~ and a + light polarization, 
respectively. Free-carrier transitions are indicated by full symbols. The filling factor was evaluated assuming somewhat higher 
2D hole density of 8.4X10 11 cm -2 as is discussed in the text. 



near 6 T. The n — level should be occupied by holes 
for filling factors v ^ 6, i.e. up to B = 5.2 T as calcu- 
lated using the nominal 2D holes density. However the 
H4 line vanishes only at B 5.8 T (see Fig. ^|). This can 
be explained assuming that the actual hole concentration 
Ns ~ 8.4 x 10 11 cm ~ 2 exceeds the nominal value as an 
effect of sample illumination. 

According to Fig. ||the line H3 observed in the a + po- 
larization corresponds to the n = — 1 hole Landau level. 
In the field range 2 T ^ B ^ 9 T this is the fifth level 
(counting from the HH1 subband edge) and thus the H3 
emission should disappear for filling factor v = 4 (i.e. 
B = 8.7 T). Indeed the peak intensity starts to fall at 
B«8T but for B > 9 T grows again and the line per- 
sists up to much higher fields B > 11 T. Furthermore 
as we can see m Fig. about B w 8.2 T the broadening 
of H3 line temporarily increases and even an additional 
peak emerges on the low-energy shoulder. This can eas- 
ily be explained as a result of the crossing of n — — 1 and 
n = 3 levels predicted by the calculations (Fig. |5|). It was 
shown (2f| that cubic terms in Hamiltonian (g) neglected 
in the axial approximation lead to strong Landau-level 
repulsion and to anticrossing behavior if the numbers n 
of levels differ by 4. This probably accounts for a slight 
discrepancy between calculated and measured energies of 
the n = — 1 level near the crossing noticeable in Fig. 

In magnetic fields where the H3 transition disappears 



a new line labeled as E2 emerges in the a + emission and 
rapidly gains in intensity. Its energy is greater than that 
of line El by the value /J-sgB of electron spin splitting. 
The detailed PL spectra are shown in Fig. |. The en- 
ergy of H3 transition can not be determined accurately 
for B > 11 T because it merges with the acceptor-bound 
exciton line A°X. However we can find the field value 
at which the n = — 1 hole level becomes depopulated by 
examining the donor-related lines D2 and D3. They cor- 
respond to recombinations of donor-bound electrons with 
n = —2 and n — — 1 holes, respectively. The line D3 dis- 
appears at B 11.6 T, which agrees well with the value 
v = 3 calculated using the hole density Ns — 8.4 x 10 11 
cm~ 2 . At higher magnetic fields the peak H2 can not be 
separated from the strong E2 line. For this reason the 
H2 energies are not plot in Fig. M. However the donor 
line D2 related to the same hole level n = — 2, remains 
visible up to the highest magnetic fields. 

Measured energies of PL lines (in particular of low- 
energy transitions HI and H2) oscillate with the mag- 
netic field. Similar behavior was frequently observed in 
n-doped structures |32, 43 , |4| [l5| and explained theo- 
retically as a result of many-body interactions fl6|| . As in 
Ref. |32| we find that the amplitude of oscillations grows 
with B and the maximal and minimal values of transi- 
tion energies correspond to even and odd filling factors, 
respectively. However in our p-doped structures the in- 
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FIG. 6: Detailed PL spectra obtained in low magnetic fields 
in the o~ polarization. Near B — 5.6 T the H4 line involving 
7i = holes declines and is replaced by a stronger peak with 
lower energy denoted as El. At the same fields a high-energy 
donor transition (presumably related to n = holes) also 
disappears. 
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FIG. 7: Parts of PL spectra in the a + polarization including 
free- hole lines measured near the crossing of n = — 1 and n — 3 
valence-band levels. Due to wave-function mixing holes from 
the level n — 3 can recombine with n — electrons and a 
weak peak (indicated by an arrow) appears on the low-energy 
shoulder of H3 transition. 
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FIG. 8: High-field PL spectra obtained in the a + polarization. 
The lines H3 and D3 related to the n = — 1 hole Landau 
level disappear at B = 11.6 T and at the same field a new 
peak named E2 becomes visible. In still higher fields the new 
transition gains in intensity and masks the line H2 related to 
n = — 2 holes. However the donor transition D2 also involving 
n = — 2 holes remains visible. 



tensity of luminescence increases at even filling factors, 
in contrast to what is reported in Ref. |32j . This can be 
seen in Fig. || where we plot the lineshape of HI transi- 
tion in the fields range 9 T < B < 14 T. The maximum 
of transition energy (y = 2) corresponds to B = 15.7 
T above the limit of fields attainable in our experiment. 
To our knowledge the origin of magnetic-field induced 
oscillations of photoluminescence intensity remains unre- 
solved. 

Apart from the oscillations the agreement between cal- 
culated and measured energies of hole Landau levels is 
very good. Also the intensities of observed transitions are 
consistent with matrix elements shown in Table II. The 
strongest transition H2 has intensity on average 2.5 — 3 
times greater than the line HI. Furthermore evaluated 
intensities of light-hole type transitions (H3, H4) are ap- 
proximately three times less than those of heavy-hole 
type ones (HI and H2). 

The origin of lines El and E2 is not clear. They rep- 



11 



c 

I! 



c5 



w 
c 




1.510 



1.515 
Energy (eV) 



1.520 



FIG. 9: Fragments of high-field PL spectra measured in the 
<j~ polarization. As the free-carrier line HI shifts to higher 
energies its intensity increases. 



resent in different circular polarizations the same opti- 
cal transition not related to a free-hole recombination. 
It is striking that the transition energy as plotted in 
Fig. H coincides with the zero-field position Ep of the 
hole Fermi level. This suggests a Fermi-edge singularity. 
However due to small masses of minority carriers in our 
structure (electrons in the GaAs layer) a possible sin- 
gularity should be quite broad (smeared over the range 
(mh/meYEp p7|). We see sharp peaks instead (Figs. 
|| and |9j) of a width below 1 meV. Moreover, as was 
observed in n-doped structures pS| , in strong magnetic 
fields the Fermi-edge singularity appears not as a new 
emission line but as some increase of the Landau-level 
transition when the level energy approaches the Fermi 
level. The other possible explanation - a recombination 
of electrons with acceptor-bound holes - must also be ex- 
cluded. The line related to this transition can be noticed 
in the spectra (see Fig. ^) as a weak feature at energy 
about 15 meV below the E-lines. 

In our view narrow widths and high intensities of the 
E-lines suggest their cxcitonic nature. Similar line was 
seen by Kim et al j32| in the PL emission from a n- 
doped single GaAs/GaAlAs heterojunction. The peak 
red-shifted with respect to the free-electron transition 
which emerged at v = 2 was also interpreted as an ex- 
citonic recombination. In our structures photoexcited 
electrons are removed by the interface electric field far 
from the junction. This effect as well as screening of 
Coulomb interaction by 2D holes prevents a formation of 
excitonic states. This explains an absence of exciton lines 
in the PL emission at low magnetic fields. Nevertheless 
they can emerge at higher fields as the strong magnetic 
field applied in the growth direction brings an electron 
closer to a hole in the interface plane and increases ex- 
citonic binding. The resonance condition between states 



at the Fermi level and those of the excitonic level pos- 
sibly further increases the intensity of E-line emission 
p9| , p0| [5lj |. It is still not clear why the excitonic peak 
appears in much higher magnetic fields in the er + than in 
the a" polarization. As we have already noticed the E- 
line emerges only when a free-hole transition disappears. 
One can try to answer this question using the concept of 
phase-space filling. According to the theory |52[ a hole 
Landau level can contribute to the excitonic wave func- 
tion as soon as it becomes depopulated by carriers. Both 
levels n — —2 and n = — 1 in the HH1 subband optically 
active in the a + polarization are occupied by holes down 
to the filling factor v — 3. At higher magnetic fields the 
empty n = — 1 states add to the excitonic level which 
can then produce also the er + radiation. We intend fur- 
ther investigations, both experimental and theoretical, to 
clarify this effect. 

In strongly doped structures with high densities of 
free 2D carriers many-body effects beyond the Hartree 
approximation affect subband positions. It was shown 
theoretically Q that in n-doped hetero junctions cor- 
rections to electron energies due to exchange-correlation 
effects are quite small. However they should be more sig- 
nificant in p-type systems because of large hole effective 
masses (in the local-density approximation the exchange- 
correlation potential is proportional to the effective Ryd- 
berg J39|). In previous calculations of hole Landau levels 
the many-body corrections were neglected. Their theo- 
retical description is difficult due to complex structure 
of the valence band edge. A solution of this problem by 
means of the local-density approximation was proposed 
by Bobbert et al. |^6| and then successfully compared 
with the results of zero-field PL measurements on Be 6- 
doped GaAs/Gai-^Al^As quantum wells We ex- 

amined the influence of many-body effects on Landau 
level energies in our structures by adding the exchange- 
correlation potential from Ref. to the Hartree poten- 
tial V(z) during self-consistent calculations. The results 
A comparison with Fig. ^ reveals 



are shown in Fig. 10 



that exchange-correlation effects change significantly the 
subband edges. The ground HH1 subband is shifted by 
12 meV to higher energies and as a result the effective 
band gap is reduced. This effect is known as a band gap 
renormalization. Furthermore the splitting between the 
HH1 and LH1 subbands is increased by about 3 meV. 
To compare calculated Landau levels with experimental 
results we once more subtracted the energies of measured 
transitions from the energy E cv + 0.5Hio c ±0.5fj,BgB, this 
time using the value E cv 12 meV greater. As is seen the 
agreement remains very good for the highest levels n = 1 
and n = —2 but is substantially worse for the others. We 
verified that the fit can not be improved by small changes 
of the hole density value used in calculations. However 
the new value of energy E cv is more consistent with ex- 
perimental results. An inspection of the conduction band 
profile in the GaAs layer (computed by adding the band 
gap energy E g — 1519.5 meV to the calculated poten- 
tial V(z) shown in Fig. ITof) reveals that it has energy 
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FIG. 10: Same as Fig. H , with exchange-correlation corrections included in the calculations. 



1449.5 meV in the flat-band region far from the inter- 
face. Hence the value E cv — 1464 meV corresponds to 
an electron placed at least 75 nm away from the junction, 
too far to recombine with confined holes. Increasing the 
energy E cv by 12 meV one shifts the electron to a rea- 
sonable distance of about 40 nm |33|. We thus conclude 
that the model of Bobbert et al. pq] adequately de- 
scribes exchange-correlation corrections to the subband 
edge but distorts the effective masses of holes in excited 
(low-energy) Landau levels. 



V. SUMMARY 

We investigated photoluminescence from p-doped sin- 
gle heterostructures at helium temperatures and in 
a magnetic field directed normally to the interfaces. 
Very low excitation powers applied in the measure- 
ments allowed to obtain rich spectra containing both free 
and bound-carrier transitions involving two-dimensional 
holes. The free-carrier lines were used to derive the struc- 
ture of valence-band Landau levels in confined subbands. 
We calculated energies and wave functions of hole states 
using a new numerical procedure based on a realistic 
model of charge distribution in studied heterojunctions. 
The theoretical results were consistent with experimen- 
tal data. All allowed optical transitions between holes in 
occupied valence-band Landau levels and photocreated 
electrons from the ground conduction-band Landau level 
were detected. Their relative intensities were coherent 
with evaluated selection rules. The measured recombi- 



nation energies agreed well with computed positions of 
hole levels. Also the field ranges in which individual 
transitions remained visible matched the calculated se- 
quence of levels. Observation of Landau-level crossing 
provided information about the hole states which - due 
to symmetry constrains - cannot normally participate to 
the emission. Detailed comparison of experimental data 
with the theory allowed discerning the optical transitions 
(named as E) with energies in the range of free-carrier re- 
combinations but not matching the Landau-level scheme. 
Interestingly the E-line positions correlated with the hole 
Fermi energy in zero magnetic field. Moreover the line 
emerged in the a + polarization in much higher magnetic 
fields than in the <j~ one. We presented some arguments 
for the excitonic nature of related optical transition pos- 
sibly influenced by a Fermi-edge singularity. Many-body 
interactions affected the PL spectra in different ways. 
We observed well-known oscillations of transition ener- 
gies in a magnetic field correlated with integer values of 
the hole filling factor. Inclusion of exchange-correlation 
effects in Landau-level calculations along with the model 
recently proposed by Bobbert et al. |26| gave mixed re- 
sults. Computed correction to the ground-subband edge 
improved the agreement with experimental results but 
we obtained incorrect values of hole masses, particularly 
in excited Landau levels. 
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